Inclusive cross section and single transverse spin asymmetry for very forward neutron 
production in polarized p+p collisions at y/s = 200 GeV 
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The energy dependence of the single-transverse-spin asymmetry, An, and the cross section for 
neutron production at very forward angles were measured in the PHENIX experiment at RHIC for 
polarized p+p collisions at ^=200 GeV. The neutrons were observed in forward detectors covering 
an angular range of up to 2.2 mrad. We report results for neutrons with momentum fraction of 
xp=0A5 to 1.0. The energy dependence of the measured cross sections were consistent with xp 
scaling, compared to measurements by an ISR experiment which measured neutron production in 
unpolarized p+p collisions at y / s=30.6-62.7 GeV. The cross sections for large xf neutron production 
for p+p collisions, as well as those in e + p collisions measured at HERA, are described by a 
pion exchange mechanism. The observed forward neutron asymmetries were large, reaching An — 
—0.08 ± 0.02 for xf—0.8; the measured backward asymmetries, for negative xf, were consistent 
with zero. The observed asymmetry for forward neutron production is discussed within the pion 
exchange framework, with interference between the spin-flip amplitude due to the pion exchange 
and nonflip amplitudes from all Reggeon exchanges. Within the pion exchange description, the 
measured neutron asymmetry is sensitive to the contribution of other Reggeon exchanges even for 
small amplitudes. 

PACS numbers: 13.85.Ni,13.88.+e,14.20.Dh,25.75.Dw 



I. INTRODUCTION 

With the first polarized p+p collisions at ^/s = 200 
GeV at the Relativistic Heavy Ion Collider (RHIC), a 
large single transverse spin asymmetry (An) for neutron 
production in very forward kinematics was discovered by 
a polarimeter development experiment [![. That exper- 
iment was designed to measure the asymmetry for very 
forward photons and used an electromagnetic calorime- 
ter. The calorimeter was used to identify neutrons, orig- 
inally to remove them from the photon data, when a 
large asymmetry was observed in forward neutrons. The 
neutron energy resolution was coarse, so no cross sec- 
tion measurement was reported. The discovery inspired 
the PHENIX experiment to use existing very forward 
hadronic calorimeters, with additional shower maximum 
detectors, to measure the neutron transverse asymme- 
try at the PHENIX interaction point at RHIC with a 
significantly better neutron energy resolution. Here we 
report the first measurements of very forward inclusive 
and semi-inclusive neutron production cross sections at 
^fs = 200 GeV and measurements of An for forward and 
backward production with improved neutron energy res- 
olution. The An is a left right asymmetry written as: 



A N = 



da^ — da^ 
da^ + dcr^ 



(1) 



for yields observed to the left when facing along the polar- 
ized proton's momentum vector, where da^ (da^) is the 
production cross section when the protons are polarized 



up (down). The An with cross section measurements for 
higher energy p+p collisions provide qualitatively new 
information toward an understanding of the production 
mechanism. 



a 




(m P ,E P ) f 



N(mn,En) 



FIG. 1: A schematic diagram of neutron production, pa — > 
nX, for the Reggeon exchange model shown with Lorentz 
invariant variables s',Q 2 and t. "a" is a proton or positron 
for p+p or e + p reactions. R indicates a Regge trajectory with 
isospin odd such as it, p, a-z and Pomeron-7r in the Regge 
theory. For pion exchange, R — ir. 



Cross sections of inclusive neutron production in un- 
polarized p+p collisions were measured at the ISR from 
y/s = 30.6 to 62.7 GeV [|||. These cross sections have 
been described using One Pion Exchange (OPE) models 
In OPE, the incoming proton emits a pion which 
scatters on the other proton as shown in Figure [TJ Kine- 
matics of the neutron are characterized by two variables, 
x f and pt defined by, 



'Deceased 

^PHENIX Spokesperson: jacak@skipper.physics.sunysb.edu 



XF = PL/PL(max) = E n COs8 n /E p ~ E n /E pi (2) 

Pt = E n sm6 n ~ x F E p n . (3) 
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where pl is the momentum component of the neutron in 
the proton-beam direction, E n and E p are energies of the 
neutron and the proton beam, and 9 n is the polar angle of 
the neutron from the beam direction which is very small 
(~mrad) for forward neutron production. The measured 
cross section showed a peak around xf ~ 0.8 and was 
found to have almost no </s dependence. OPE models 
gave a reasonable description of the data. 

OPE models were also used to describe proton and 
photon induced production of neutrons measured at the 
HERA e+p collider [Til. Il2j]. These measurements probe 
the pion structure function at small x. The NA49 col- 
laboration also published the cross section for forward 
neutron production for p+p collisions at y/s = 17.3 GeV 
[HI . They compared the result with those from ISR and 
HERA and found they did not agree. 

The neutron asymmetry provides a new tool to probe 
the production mechanism. For the OPE model, An 
arises from an interference between spin-flip and spin- 
nonflip amplitudes. Since the pion-exchange amplitude 
is fully spin-flip, the asymmetry is sensitive to other 
Reggeon exchange amplitudes which are spin-nonflip, 
even for small amplitudes. 

This paper presents the xf dependence of cross sec- 
tions, inclusive and semi-inclusive (with a beam-beam 
interaction requirement), and An for very forward and 
very backward neutron production in polarized p+p col- 
lisions at = 200 GeV. 



II. EXPERIMENTAL SETUP 

A. Detector apparatus 

A plan view of the experimental setup for very forward 
neutron measurement at PHENIX [3] is shown in Fig. [2] 
The RHIC polarized proton beams were vertically polar- 
ized. Each collider ring of RHIC was filled with up to 
111 bunches in a 120 bunch pattern, spaced 106 ns apart, 
with predetermined patterns of polarization signs for the 
bunches. The colliding beam rotating clockwise when 
viewed from above is referred to as the "Blue beam" and 
the beam rotating counterclockwise, the "Yellow beam" . 
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Dx 



Dx 
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Collision point 



Yellow Beam ~ Blue Beam 

□ ZDC (W-Cu alloy) □ Charge veto counter (Plastic Scintillator) 
13 SMD (Plastic Scintillator) 

FIG. 2: A plan view of the experimental setup at PHENIX 
(not to scale). Shown are the principal components for the 
neutron physics. Charged veto counters are in front of ZDCs, 
and the SMDs are between the first and second ZDC modules. 

Neutrons were measured by a Zero-Degree Calorimeter 
(ZDC) [lH with a position-sensitive Shower-Maximum 



Detector (SMD). One ZDC module is composed of Cu-W 
alloy absorbers with PMMA-based communication grade 
optical fibers, and corresponds to 1.7 nuclear interaction 
lengths. A single photomultiplier collects Cerenkov light 
via optical fibers. Three ZDC modules are located in 
series (5.1 nuclear interaction lengths) at ±1800 cm away 
from the collision point, covering 10 cm in the transverse 
plane. 

The SMD comprises x-y scintillator strip hodoscopes 
and is inserted between the first and second ZDC modules 
(see Fig. 5 of [Hj]) at approximately the depth of the max- 
imum of the hadronic shower. The x-coordinate (horizon- 
tal) is given by 7 scintillator strips of 15 mm width, while 
the y-coordinate (vertical) is given by 8 strips of 20 mm 
width, tilted by 45 degrees. 

The neutron position can be reconstructed from the en- 
ergy deposited in scintillators with the centroid method. 
We calculated the centroid: 



E(i) ■ x(i) 



ArSMD 



(4) 



where E(i) and x(i) are the energy deposit and the posi- 
tion of the i-th scintillator, respectively. The number of 
scintillators with pulse height above the Minimum Ion- 
ization Particle (MIP) peak is shown as which is 
defined as the SMD multiplicity. 

Detectors are located downstream of the RHIC dipole 
(DX) magnet, so that collision-related charged particles 
are swept out. A forward scintillation counter, with di- 
mensions 10 x 12 cm, was installed in front of the ZDC to 
remove charged particle backgrounds from other sources. 
In this analysis, we used only the south ZDC detector, 
which is facing the Yellow beam. 

As a beam luminosity monitor, Beam Beam Counters 
(BBCs) are used. The BBC comprises 64 photomultiplier 
tubes and 3 cm thick quartz Cerenkov radiators. The 
two BBCs are mounted around the beam pipe ±144 cm 
away from the collision point which cover ±(3.0-3.9) in 
pseudorapidity and 27r in azimuth. 

The neutron data were collected in 2006 with two trig- 
gers. One is the ZDC trigger for neutron inclusive mea- 
surements, requiring an energy deposit in the south ZDC 
greater than 5 GeV. The other trigger was a ZDC<K>BBC 
trigger, a coincidence trigger of the ZDC trigger with 
BBC hits which are defined as one or more charged par- 
ticles in both of the BBC detectors. We note that the 
ZDC trigger was prescaled due to data acquisition limi- 
tations. Therefore, the ZDC trigger samples are signifi- 
cantly smaller than the ZDC<X>BBC trigger samples. 



B. Detector performance 

In order to evaluate the detector performance, sim- 
ulation studies were performed with GEANT3 with 
GHEISHA [l6j] which simulated the response of the pro- 
totype ZDC to hadrons well. A single neutron event 
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generator and PYTHIA (version 6.220) [T7| were used to 
generate events. The single neutron event generator sim- 
ulated neutrons as a function oi xf and pr- The xf 
distribution which was used for the simulation input was 
determined as a differential cross section, da/dxF, in the 
cross section analysis (section Till A|) . The pt distribution 
is difficult to determine by the PHENIX data alone since 
the position and energy resolutions are insufficient to ade- 
quately determine it, so the px distribution from the ISR 
result, exp(— 4.8 px (GeV/c)), was used as simulation in- 
put, assuming pt scaling from the ISR to the PHENIX 
energies. To check the reliability of this assumption, dis- 
tributions of radial distance from the detector center, r, 
for the data and simulation were compared based on the 
relation of pt oc r as, 



p T = E n six\9 n = E n 



yV2 



En 



(5) 



where d is the distance from the collision point to the 
detector, corresponding to 1800 cm, and r is determined 
for the shower centroid with Eq. (|4|. 

The comparison of r distributions with the integration 
of measured ZDC energies 20-120 GeV agreed well as 
shown in Fig. [3] 
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FIG. 3: (color online) r distributions for the data and sim- 
ulation with the exponential pt shape. Distributions agreed 
within r<4 cm. 



1. Performance of the energy measurement 

The neutron energy measurement with the ZDC was 
degraded by a nonlinearity of the photoelectron yield and 
shower leakage out the back and sides of the detector 
(edge effect). The ZDC response was studied by simula- 
tion with the single neutron event generator. 

The energy linearity and resolution were evaluated 
from the response to incident neutrons with energies from 



20 to 100 GeV in the simulation. The absolute scale was 
normalized at 100 GeV with the experimental data. The 
ZDC response below 100 GeV exhibits nonlinear behavior 
as shown in Fig. @] We applied a correction of the non- 
linearity to the experimental data based on this result. 
We used the difference between the linear and nonlinear 
response as a component of the systematic uncertainty 
in the determination of the cross section (section IHI A[) . 
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FIG. 4: (color online) Top) The mean of output energy as 
a function of the incident neutron energy evaluated by the 
simulation. Solid line indicates a linear response. Bottom) 
The energy resolution as a function of l/y/E (GeV" 1/2 ). Solid 
line shows the fit result; AE/E = 65%/VE + 15%. 

As shown in Fig. @J the energy resolution for 20-100 
GeV neutrons was described by 



AE 
~E~ 



65% 



VEjGeV) 



15%. 



(6) 



The absolute scale of the energy measurement was nor- 
malized with the 100 GeV single neutron peak in heavy 
ion collisions. However, the energy of neutrons from p+p 
collisions was below 100 GeV, so simulation was used to 
estimate the detector response for neutron energies in 
this region. 

Figure [5] shows the absolute energy scale calibrated 
by observing one neutron from peripheral Cu+Cu colli- 
sions at y/SNN = 200 GeV; 100 GeV neutrons less than 
2 mrad from the beam axis produced the single neutron 
peak. The energy resolution expected from simulation 
was about 22% for the 100 GeV neutron and was consis- 
tent with the observed width of the single neutron peak 
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as shown in Fig. [5j The energy nonlinearity was con- 
firmed by the single neutron peak from Cu+Cu collisions 
at y/SNN = 62.4 GeV shown in Fig. [5] which peaked at 
26±0.3 GeV, consistent with nonlinearity indicated by 
the simulation. 



prototype reproduced this effect. 

A residual edge effect was seen in the data at the top 
and bottom of the detector, so we chose to apply a fidu- 
cial cut to minimize the effect. According to the sim- 
ulation, 95-100% of the incident energy was contained 
within r < 3 cm. 
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FIG. 5: The energy distribution in the ZDC for Cu+Cu col- 
lisions at ^snn = 200 GeV. Peripheral events were selected 
by requiring BBC inactivity. 



2. Performance of the position measurement 

The position resolutions were evaluated by the simu- 
lation. Figure [7] shows the position resolution (RMS) as 
a function of the neutron incident energy for x (horizon- 
tal) and y (vertical) positions. The position resolution 
was approximately 1 cm for the neutron energy at 100 
GeV. 
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FIG. 6: (color online) The energy distribution in the ZDC 
for Cu+Cu collision at ^/snn = 62.4 GeV. The neutron peak 
position was determined with a Gaussian + polynomial fit. 

The edge effect was studied by a prototype ZDC with a 
100 GeV proton beam at CERN. Generally, the measured 
energy decreased near the edge, however, nearest the 
PMT, the measured energy increased. This was found 
to be caused by the fibers in the top region which con- 
nected to the PMT (see Fig. 5 of [la]); where the shower 
hit the fibers directly. The simulation used to study the 
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FIG. 7: (color online) The position resolution (RMS) as a 
function of the incident neutron energy for x (top) and y 
(bottom). Circles show measured values. They were well 
reproduced by a second order polynomial fit. Red lines show 
the fit results. 

Near the edge of the detector, the position measure- 
ment is also affected by shower leakage. If the incident 
position was in the edge area, the output position was 
shifted to the detector center due to shower leakage, in- 
dependent of neutron energy. This position shift caused 
by the edge effect is corrected based on the simulation. 

The reliability of the position measurement was stud- 
ied by comparing hadron shower shapes of the data and 
simulation. The shower width and highest shower frac- 
tion among all scintillators were calculated for x and 
y independently. We compared the measured distribu- 
tion with simulation for each SMD multiplicity since the 
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hadron shower shape sensitively depends on the SMD 
multiplicity. The distribution of y was well reproduced by 
the simulation, however the distribution of x was not well 
reproduced, especially for the highest shower fraction in 
high SMD multiplicity events. The systematic uncer- 
tainties for the position measurement were estimated by 
matching the highest shower fraction of x by smearing the 
simulated shower shapes in case of the SMD multiplicity 
= 7, which shows the worst agreement between the data 
and simulation. After the smearing to match the highest 
shower fraction, the shower width of the simulation also 
reproduced that of the experimental data. The position 
resolution increased 14% after the smearing. 



3. Performance of the neutron identification 

Events within the detector acceptance in p+p collisions 
were studied with GEANT3 with PYTHIA event generators, 
and the performance of neutron identification and its re- 
liability were evaluated. 

We studied particle species detected in the ZDC with 
the 5 GeV energy threshold which was required for the 
ZDC trigger (without the BBC coincidence requirement). 
In about 92% of events, only a single particle was de- 
tected by the ZDC in each p+p collision, mainly pho- 
tons, neutrons and protons. Energy distributions for 
these three particles are plotted in Fig. [5J 
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FIG. 8: (color online) Energy distributions in ZDC for neu- 
tron, photon and proton. The ZDC threshold was set at 5 
GeV in the ZDC trigger. Events with one particle detected 
in the ZDC in each p+p collision are shown. 

Only neutral particles, photons and neutrons, were ex- 
pected to be detected with the ZDC due to sweeping of 
charged particles by the DX magnet. However scattered 
protons could hit the DX magnet or beam pipe and cre- 
ate a hadronic shower and particles from the shower could 
hit the ZDC. 



Most of the photons and neutrons were generated by 
diffractive and gluon scattering processes. In pythia 
hard processes, neutrons are generated mainly from 
string fragmentation (~65%) and then decay from 
A°,A + ,A~,A°. The forward photons were generated 
by decays of ir°s (~91%) and ijs (~7%). Protons were 
generated by elastic and diffractive processes. Particles 
depositing less than 20 GeV of energy in the ZDC were 
predominantly photons and protons as shown in Fig. [8j 

Photons are mostly absorbed in the first ZDC mod- 
ule, which is 51 radiation lengths long. Thus, photons 
were removed by requiring energy deposited in the SMD 
or in the second ZDC module. In photon rejection with 
the SMD, more than one scintillator above threshold (the 
SMD multiplicity > 2) were required for both x and y. 
After applying this cut, the neutron purity was estimated 
to be 93.6±0.3%. In photon rejection with the second 
ZDC module, energy deposited in the second ZDC mod- 
ule above 20 GeV was required. After applying this cut, 
the neutron purity was estimated to be 93.6±0.5%. In 
the analyses of the cross section and the asymmetry, pho- 
ton rejection with the SMD was applied since the position 
information calculated by the SMD was required. Rejec- 
tion with the second ZDC was used for the estimation of 
the rejection efficiency with the SMD which is discussed 
in section UlI Al 

The charge veto counter was used to reject protons. A 
neutron energy above 20 GeV and the charge veto cut 
removed most proton events, as discussed later in this 
section. 

The main backgrounds after neutron identification are 
K°s and protons. The purities were estimated for neu- 
tron energies above 20 GeV. In the cross section and the 
asymmetry analyses, we also required the acceptance cut 
and/or a higher energy cut. In these cases, the purities 
improve and are estimated in each analysis section. 

In the ISR experiment, the K contamination to the 
neutron measurement was estimated from the K mea- 
surements ||. They obtained 10 % contamination at 
Xi?=0.2 and less than 4 % at xf >0.4. The fraction 
of K° to neutron in pythia is consistent with the ISR 
result. We have included no correction for the K° con- 
tamination in this analysis. 

The proton background is very sensitive to the ma- 
terials around the ZDC and the magnet tuning in the 
accelerator. The systematic uncertainty of proton con- 
tamination was estimated by the simulation using the 
measured fraction of charged events in the charge veto 
counter. Noise was estimated by the pedestal width of 
the data and was incorporated into the simulation. For 
the proton contamination analysis, photon events were 
removed by requiring the second ZDC module cut. The 
fraction of proton events can be estimated as a fraction of 
charged candidates, which are the events with one more 
MIPs in the charge veto counter. These fractions were 
0.42 and 0.28 for the data and simulation, respectively. 
Proton events in the experimental data were about 1.5 
times more frequent than that of the simulation. We 
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ascribe the difference to beam conditions that cause in- 
teractions with materiais around the DX magnet and the 
ZDC. 

The threshold dependence of the selection of charged 
particle candidates was also studied. The change in 
charged particle fraction was less than 1% so that the 
threshold dependence was negligibly small. Therefore, 
the factor 1.5 was a reasonable estimate for the fraction 
of charged candidates between the data and simulation. 

The proton background was estimated and included 
in the systematic uncertainties. According to the sim- 
ulation study for the structure of proton events, proton 
events should be detected in the direction of beam bend- 
ing which is negative x for the south ZDC. This behav- 
ior was confirmed by the experimental data as shown in 
Fig. IHl which is a plot of the x position determined by 
the SMD vs. the charge distribution in the charge veto 
counter. Most charged candidates were distributed in the 
negative x region. We assumed the proton background 
An equals zero and evaluated its systematic uncertainty 
by the dilution method with ^4^=0 in the asymmetry 
analysis. 
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FIG. 9: The x position calculated by the SMD vs. the charge 
distribution in the charge veto counter for the experimental 
data. Most charged events were distributed in the negative 
x region which is the direction of beam bending by the DX 
magnet. 



III. CROSS SECTION MEASUREMENT 

A. Analysis 

The differential cross section with respect to xf was 
measured: 

du _ -Ancutron 1 

dx f £ dxp' 

where A ncutr0 n is the number of neutrons after the cor- 
rection of cut efficiencies and the energy unfolding. 



For the cross section analysis, 6.5 million events taken 
by the ZDC trigger were used from the sampled luminos- 
ity of 240 nb _1 . The acceptance cut r<2 cm was used 
to select kinematics similar to the ISR experiment. We 
assumed the beam axis was the same as the ZDC cen- 
ter in this analysis and the deviation was evaluated as a 
systematic uncertainty. The beam axis compared to the 
ZDC center is discussed in Appendix A. 




ZDC energy (GeV) 

FIG. 10: The energy distribution measured with the ZDC 
after the neutron identification and the acceptance cut (r<2 
cm, corresponds to pt<0.11 GeV/c). 

Figure [10] shows the energy distribution measured with 
the ZDC after the neutron identification and the accep- 
tance cut. The energy spectrum was peaked at about 70 
GeV, and this was used for a stability check of the ZDC 
gain run by run, which was found to be stable. The ratio 
of the neutron yield to the BBC counts without the colli- 
sion vertex requirement was used for a stability check of 
the neutron selection, and it was also found to be stable. 

One background source was beam-gas interaction. 
Beam-gas events are normally removed by requiring 
a forward-backward coincidence of the BBC detectors. 
However, this could not be done for the ZDC triggered 
events. Instead, we evaluated the fraction of beam-gas 
background using the 9 noncolliding bunch crossings with 
the combination of filled and empty bunches at PHENIX. 
We found that the fraction was 0.0062±0.0004 on aver- 
age, negligibly small. 

The neutron hit position was calculated by the cen- 
troid method using the distribution of scintillator charge 
above the threshold in the SMD, Eq. Q. In this analysis, 
the same threshold was applied to the data and simula- 
tion and the efficiency of the SMD cut was estimated by 
simulation. The difference of efficiencies caused by uncer- 
tainty of the SMD cut efficiency was estimated using the 
nearly pure neutron sample by the neutron identification 
with the second ZDC cut (section III F3 3|) . The energy 
spectrum was corrected based on the SMD cut efficiency 
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before the energy unfolding. 

The measured neutron energy with the ZDC is smeared 
by the energy resolution. To extract the initial energy 
distribution, it is necessary to unfold the measured en- 
ergy distribution. The energy unfolding method is de- 
scribed in Appendix B. 

The ZDC energy response to neutrons below 100 GeV 
was found by the simulation to be nonlinear as described 
in section III B II This nonlinearity was included in the 
transition matrix A of Appendix B, and corrected by the 
energy unfolding. Since the hadronic interaction could 
only be determined from simulation, a systematic un- 
certainty was included, using the variation of the cross 
section evaluated with a different matrix A with a linear 
response. 

The efficiency of the experimental cuts, including the 
neutron identification and the acceptance cut, for the 
unfolded x f distribution was estimated by the simulation 
with the single neutron event generator. The acceptance 
cut used the radius, r, and the efficiency was evaluated 
from the pr distribution based on Eq. ([3]) : px ~ xp ■ E p ■ 
9 n « 0.056 -x F -r GeV/c. 

For the pr distribution, we used two distributions: a 
Gaussian form da/dpr oc exp(— ap^), where ci(xf) was 
obtained by HERA [11] with error evaluation, and an ex- 
ponential form da/dpT oc exp(— bpx), where b = 4.8±0.3 
(GeV /c) ~ 1 which was used in the ISR analysis @, Q . The 
simulated pt distributions with those two input distribu- 
tion were compared with experimental data normalized 
to the same total entries. It was found that the differ- 
ences between data and those two inputs were not large 
as shown in Fig. ITT1 
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FIG. 11: (color online) Comparison of the pr distribution 
from experimental data (black closed circles) and two simula- 
tions using Gaussian form (blue open squares) and exponen- 
tial form (red open circles) inputs. 

Figure IT21 shows the simulated pt distributions (dashed 
line) in each xf bin. The geometrical maximum px for 



the acceptance cut, r < 2 cm, in each xf are given by 
p^J ax = 0.11 • x p GeV/c, shown as dot-dashed vertical 
lines. The actual pr distributions with the experimental 
cuts were smeared due to the position resolution and the 
energy resolution, shown as solid lines. Ratios of these 
counts are the efficiency for the experimental cuts, and 
are listed in Table |TJ The errors were derived considering 
the uncertainty in the parameter a(xp) in the Gaussian 
form evaluated by HERA. There is no significant differ- 
ence in the result in case of using the ISR (exponential) 
Pt distribution. 
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FIG. 12: Simulated pr distributions using the Gaussian pr 
shape. Initial pr distributions are shown as dashed lines in 
each xf region. The expected pr region for the acceptance, 
r < 2 cm, is below the vertical dot-dashed lines, which are the 
maximum pr calculated as w 0.11 ■ xf GeV/c. The actual pr 
distributions with the experimental cuts are shown as solid 
lines. 



TABLE I: The expected pr for r <2 cm, mean pr value with 
the experimental cut, and the efficiency for the experimental 
cut estimated by the simulation (Fig. I12|) . The errors were 
derived considering the uncertainty in the parameter cl(xf) 
in the Gaussian form evaluated by HERA. 



neutron 


mean 


efficiency 


xf 


Pt 






(GeV/c) 




0.45-0.60 


0.072 


0.779 ± 0.014 (1.8%) 


0.60-0.75 


0.085 


0.750 ± 0.009 (1.2%) 


0.75-0.90 


0.096 


0.723 ± 0.006 (0.8%) 


0.90-1.00 


0.104 


0.680 ± 0.016 (2.3%) 
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The mean values of the simulated px distributions in 
each energy region are also listed in Table HI The cross 
section was obtained after the correction of the energy 
unfolding and the cut efficiency 



TABLE II: Systematic uncertainties for the cross section 
measurement. The absolute normalization error is not in- 
cluded in these errors. The absolute normalization uncer- 
tainty was estimated by BBC counts to be 9.7% (22.9T2.2 
mb for the BBC trigger cross section). 

exponential pr form Gaussian pr form 
3 - 10% 7 
3 - 31% 



Pt distribution 
beam center shift 
proton background 
multiple hit 

total 



22% 



3.6% 
7% 



11 - 33% 



16 - 39% 



Table ITll summarizes all systematic uncertainties evalu- 
ated as the ratio of the variation to the final cross section 
values. The absolute normalization error is not included 
in these errors. It was estimated by BBC counts to be 
9.7% (22.9±2.2 mb for the BBC trigger cross section). 

The background contamination in the measured neu- 
tron energy with the ZDC energy from 20 to 140 GeV 
for the acceptance cut of r < 2 cm was estimated by 
the simulation with the pythia event generator. The 
background from protons was estimated to be 2.4% in 
the simulation. The systematic uncertainty in the exper- 
imental data was determined to be 1.5 times larger than 
this as discussed in section Hi B 31 Multiple particle de- 
tection in each collision was estimated to be 7% with the 
r < 2cm cut. 

In the cross section analysis, we evaluated the beam 
center shift described in Appendix A as a systematic un- 
certainty. For the evaluation, cross sections were calcu- 
lated in the different acceptances according to the result 
of the beam center shift while requiring r<2 cm, and the 
variations were applied as a systematic uncertainty. 



Result 



TABLE III: The result of the differential cross section 
da/diF (mb) for neutron production in p+p collisions at 
V / s=200 GeV. The first uncertainty is statistical, after the 
unfolding, and the second is the systematic uncertainty. The 
absolute normalization error, 9.7%, is not included. 



(x F ) 


exponential px form 


Gaussian pr form 


0.53 


0.243±0.024±0.043 


0.194±0.021±0.037 


0.68 


0.491±0.039±0.052 


0.455±0.036±0.085 


0.83 


0.680±0.044±0.094 


0.612±0.044±0.096 


0.93 


0.334±0.035±0.111 


0.319±0.037±0.123 
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FIG. 13: (color online) The cross section results for forward 
neutron production in p+p collisions at y / s=200 GeV are 
shown. Two different forms, exponential (squares) and Gaus- 
sian (circles), were used for the pr distribution. Statistical 
uncertainties are shown as error bars for each point, and sys- 
tematic uncertainties are shown as brackets. The integrated 
Pt region for each bin is < pr < O.Uxf GeV/c. Shapes of 
ISR results are also shown. Absolute normalization errors for 
the PHENIX and ISR are 9.7% and 20%, respectively. 

The differential cross section, da/dxF, for forward neu- 
tron production in p+p collisions at y/s=200 GeV was 
determined using two pt distributions: a Gaussian form, 
as used in HERA analysis, and an exponential form, used 
for ISR data analysis. The results are listed in Table IIIII 
and plotted in Fig. [13] We show the results for xp above 
0.45 since the data below 0.45 are significantly affected by 
the energy cut-off before the unfolding. The px range in 
each xf bin is < pt < O.Uxf GeV/c from Eq. @ with 
the acceptance cut of r < 2 cm. The absolute normaliza- 
tion uncertainty for the PHENIX measurement, 9.7%, is 
not included. 

Invariant cross sections measured at the ISR exper- 
iment were converted to differential cross sections for 
the comparison with the PHENIX data. The conver- 
sion formula from the invariant cross section Ed 3 a /dp 3 
to da/dxF is described with the approximation in the 
forward kinematics as 



da 2tt f d 3 a 
- — = — / E—p T dp T , 

dx F X F Jacc d P 



(8) 



where Acc. means the pp range of the PHENIX accep- 
tance cut; < pt < O.Uxf GeV/c for the r < 2 cm cut. 
As apx shape, we used an exponential form exp(— A.8pp) 
which was obtained from the 0.3 < xp < 0.7 region from 
the ISR results @,[|. 
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For both the table and figure, we give the PHENIX 
results for two pt shapes, the exponential shape used for 
the ISR results, and the Gaussian shape used for HERA 
results. 

The measured cross section at y/s—200 GeV is consis- 
tent with the ISR result, indicating that xp scaling is 
satisfied at the higher center of mass energy. This result 
is consistent with the OPE model. 



IV. SINGLE TRANSVERSE SPIN 
ASYMMETRY MEASUREMENT 

A. Analysis 

The single transverse spin asymmetry is obtained from 
the azimuthal modulation of neutron production relative 
to the polarization direction of a transversely polarized 
beam on an unpolarized target, and normalized by an in- 
dependent measurement of the beam polarization. The 
stable polarization direction of protons is vertical with 
respect to the accelerator plane. There is an approx- 
imately equal number of bunches filled with the spin 
of polarization-up protons as of polarization-down pro- 
tons. With both beams polarized, single-spin analyses 
were performed by taking into account the polarization 
states of one beam, averaging over those of the other. 
The beam polarizations were measured using fast carbon 
target polarimeters [l8[ at a different location at RHIC 
with several measurements in each fill. The carbon tar- 
get measurements were normalized to absolute polariza- 
tion measurements made by a separate polarized atomic 
hydrogen jet polarimeter |19j . The polarizations ranged 
from 0.43 to 0.48 for the Blue beam and from 0.46 to 
0.52 for the Yellow beam. Systematic uncertainty for the 
Blue beam polarization is 5.9%, and that for the Yellow 
beam polarization is 6.2%. 

The acceptance definition for the azimuthal angle (</>) 
of neutron production is shown in Fig. Qj] where the 
polarization-up direction points to (f> = 0. The accep- 
tance cut at the ZDC required 0.5<r<4.0 cm. The accep- 
tance area was divided into 16 slices in a radial pattern. 
For the asymmetry calculation, we used a square-root for- 
mula which cancels many systematic uncertainties, such 
as detector and luminosity asymmetries: 



(9) 



where iVj (N^) is the number of events with polarization- 
up (-down) producing neutrons to azimuthal angle <f>. 

A correction is applied, discussed later, to account 
for smearing from position resolution. After normaliza- 
tion by the polarization, P, we obtain the asymmetry 
as, 
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FIG. 14: The acceptance definition for the cj> dependence 
of ejv, shown as a plot of the measured neutron position at 
the ZDC. The acceptance was divided into 16 slices in a radial 
pattern and the asymmetry was calculated by the square root 
formula starting at <j> = —tv/2 to n/2. 



For this analysis, we used 6.5 million and 17.6 million 
events for the ZDC trigger sample and ZDC<8>BBC trigger 
sample respectively from the sampled luminosity of 240 
nb _1 . A ZDC energy cut was required to select 40-120 
GeV in the measured energy. 

The raw measured asymmetry ejv(</>) divided by the 
polarization are fitted to a sine: 



A((j>) = A N sin(> - O ), 



(11) 



M4>) = ^ejv(# 



(10) 



where <po allows a deviation of the maximum asymmetry 
axis from vertical. 

In the present analysis, we used only the south ZDC 
detector, which faces the Yellow beam. The forward 
neutron asymmetry uses the polarized Yellow beam and 
sums over the polarization states of the Blue beam 
bunches. Following the Basel (Ann Arbor) Convention 
[20| . a positive An indicates more production to the left 
of the polarized (Yellow) beam, for the polarization-up 
bunches in the Yellow beam. The asymmetry for neu- 
trons produced backward was measured using a polarized 
Blue beam, summing over the polarization states of the 
Yellow beam bunches. In order to follow the Basel Con- 
vention, signs of the backward An were inverted from the 
fitting results. A positive An would indicate more neu- 
tron production to the left of the Blue (polarized) beam 
for polarization-up bunches. 

We performed two sets of simulations to estimate the 
smearing parameters, C^, which were correlated to the 
neutron energy-dependent position resolution (section 
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III B ip . The energy distributions for the simulation in- 
puts were determined in the same way as the cross section 
analysis (Section fill Ap . 

The ejv(0) was smeared from the A{<f) due to position 
resolution. From Eq. (jlOp , the smearing parameter, C$, 
can be evaluated from simulation as, 

^Output 

C * = A Input ' ( 12 ) 
iV 

where A N utput corresponds to the £at(0) of the experi- 
mental data; it includes effects of the experimental cut 
and the position resolution. As J 4^ lp " t , we generated 
neutrons with the sine modulated A(4>) as Eq. (fTTj) with 
Apf=A I j^" put =— 0.10. The smeared amplitude was ob- 
tained as A N utput — —0.076 and their ratio, 0.76, is the 
correction factor of the smearing effect, — 0.760 ± 
0.015 (ZDC trigger). For the ZDC«)BBC trigger we 
obtained the smearing parameter = 0.746 ± 0.016 
(ZDC® BBC trigger). 

For the analysis of the xf dependence of An, we chose 
bins of 40-60, 60-80, and 80-120 GeV in the measured 
ZDC energy. Events with ZDC energy greater than 120 
GeV were eliminated from this analysis (3.8% of the 
events) . Similar simulations and calculations of were 
performed for the analysis of the xf dependence of the 
asymmetry with both the ZDC trigger and ZDC<g>BBC 
trigger. 

After correction for the smearing effect, we obtain the 
measured energy dependence of An- The mean xf values 
for the ZDC trigger sample and ZDC<X>BBC trigger sam- 
ple were evaluated by the simulations which were mod- 
ified to reproduce the measured energy distributions for 
each trigger sample. 

The background contamination was studied by the 
simulation with the PYTHIA event generator. In the anal- 
ysis of the xf dependence of A n, an acceptance cut of 
r < 3 cm was applied. 

After the neutron identification and the acceptance 
cut, as described in section III B 31 the neutron purities 
were 0.975 ± 0.006 for the ZDC trigger sample, and 0.977 
± 0.010 for ZDC®BBC trigger sample. Main background 
contributions were the K° and proton. According to the 
discussion in section lll B 31 we applied the systematic un- 
certainty contributed from the proton only. They were 
1.4% and 1.0% for the ZDC trigger and ZDC<g)BBC trig- 
ger respectively and were increased by the factor 1.5 es- 
timated higher frequency of proton background in the 
experimental data, compared to simulation, to give 2.1% 
and 1.5%, which were included as systematic uncertain- 
ties. Multiple particle detection in each collision was es- 
timated to be 6.5% for the ZDC trigger and 5.9% for the 
ZDC<g)BBC trigger for the 0.5<r<4.0 cm cut. 

To evaluate the systematic uncertainty for determina- 
tion of the beam axis, An were calculated with center 
positions as (x,y) = (0.46, 0.00), (0.00, -1.10) and (0.46, 
-1.10) cm while keeping the acceptance cut, 0.5<r<4.0 



cm. These values were chosen based on measurements 
of the beam center as discussed in Appendix A. Max- 
imum variations to final values, which were calculated 
by (x,y) = (0.00, 0.00) cm, were included as systematic 
uncertainties. 

Since the smearing effect was caused by the position 
resolution, the systematic uncertainty of the position res- 
olution, 14% ( section III B 2|) . should be reflected in the 
uncertainties for the result. This was evaluated with 
a variation of the asymmetry calculated with 14% in- 
creased position resolution uncertainty in the simulation. 
The asymmetry was reduced by 4.2%. This was assigned 
as a systematic uncertainty of the smearing correction. 

A technique called "bunch shuffling" was used to check 
for systematic effects in the asymmetry measurements 
due to a variation of beam characteristics bunch by 
bunch. By randomly assigning bunch polarization di- 
rections, we create data sets of experimental data with 
little or no net polarization, and compare the result- 
ing measured asymmetry with statistical uncertainties. 
The fluctuation of measured asymmetries should corre- 
spond to the statistical uncertainty. We concluded that 
the fake asymmetry from bunch characteristics is less 
than 0.39a s tat and 0.36a s t a t for the ZDC trigger and 
ZDC<g>BBC trigger respectively. We do not include these 
uncertainties in the final systematic uncertainties for An- 

PT-correlated uncertainties from the beam center shift 
were evaluated in a similar way to the cross section anal- 
ysis described in Appendix A. They were 0.004 in the Xf~ 
integrated analysis, and 0.004-0.010 in the ^-dependent 
analysis. 

TABLE IV: Scale uncertainties for the An measurements. 

ZDC trigger ZDOgiBBC trigger 

proton background 2.1% 1.5% 

multiple hit 6.5% 5.9% 

smearing 4.2% 

total 8.0% 7.4% 



Scale uncertainties are summarized in Table HVl for the 
An measurements. Values are presented as scale varia- 
tions to the final values. Total uncertainties were calcu- 
lated by quadratic sum. The scale uncertainty from the 
beam polarization is not included in the table. The un- 
certainty in the Yellow beam polarization which was used 
in the forward neutron asymmetry measurement was ± 
6.2%, and that in the Blue beam polarization used in the 
backward neutron asymmetry measurement was ± 5.9%. 

B. Azimuthal modulation of forward neutron 
production 

In this section we present the results for the azimuthal 
modulations for neutron production, within the accep- 
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tance from 6 n — 0.3 mrad (r — 0.5 cm at ZDC) to 
9 n = 2.2 mrad (r = 4 cm), and ZDC energy from 40 
GeV to 120 GeV. 
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FIG. 15: (color online) Results for the azimuthal modulation 
for forward neutron production from polarized p+p collisions 
at ^=200 GeV in the ZDC trigger sample (top) and the 
ZDGgiBBC trigger sample (bottom). The error bars indi- 
cate the statistical uncertainties. Results for a sin(</>) fit to 
the data are indicated. The pr-correlated systematic uncer- 
tainties from the beam center shift, and scale uncertainties 
listed in Table [TVl and polarization scale uncertainties are not 
included. 

Asymmetries A(4>) were calculated for eight azimuthal 
angle bins, using Eq. §§§ and Eq. (|TU)) . Figure IT51 and [TBI 
present A(<ft) for the two trigger conditions, for forward 
and backward neutron production respectively. Statis- 
tical uncertainties are shown in the figure. The pt- 
correlated systematic uncertainties from the beam center 
shift are not shown. In addition, there are scale uncer- 
tainties listed in Table IIVI and polarization scale uncer- 
tainties. 

A significant asymmetry is present for forward neu- 
tron production. The A(cf>) data were fitted with a sine 
curve, Eq. lfTTj) , to obtain ^4^. The azimuthal offsets, 
4>q, were consistent with <po = 0. The results obtained 
for A N are: A N = -0.061 ± 0.010(siai) ± O.OOA(syst) 
(% 2 /ndf = 3.05/6) for the ZDC trigger sample and An — 
-0.075 ± 0.004(siat) ± 0.004(sysi) (x 2 /ndf = 2.22/6) 
for the ZDC<g)BBC trigger sample. There is no observed 
asymmetry for backward neutron production. The re- 
sults for backward neutron production for An are: An = 
-0.006±0.011(siai)±0.004(s?/si) (x 2 /ndf = 5.18/6) for 
the ZDC trigger sample and A N = -0.008±0.005(sia£)± 



-e- 
< 




FIG. 16: (color online) Results for the azimuthal modulation 
for backward neutron production from polarized p+p colli- 
sions at ^=200 GeV in the ZDC trigger sample (top) and 
the ZDC® BBC trigger sample (bottom). The error bars in- 
dicate the statistical uncertainties. Results for a sin(0) fit to 
the data are indicated. The pr-correlated systematic uncer- 
tainties from the beam center shift, and scale uncertainties 
listed in Table [TVl and polarization scale uncertainties are not 
included. 



0.004(syst) {x 2 /ndf = 3.31/6) for the ZDOgiBBC trigger 
sample. 

To compare with the previous result [ij from the po- 
larimeter development experiment at RHIC, we com- 
pared to the An of the forward ZDC<g>BBC trigger sam- 
ple. The amplitude of the measured An was; An = 
(-0.090±0.006±0.009) x (l.Q0±^;||). Errors indicate the 
statistics, systematics and the scaling uncertainty from 
the polarization measurement. The two results are con- 
sistent within the errors, including the scaling uncertain- 
ties. We note that the two measurements used slightly 
different detection coverages for the charged particle in- 
teraction trigger: 2 .2 < 1 77 1 <3.9 in the horizontal and ver- 
tical directions for the polarimeter development experi- 
ment, and 3.0<|?7|<3.9 for the PHENIX experiment. 



C. if dependence of An 

The xf dependence of An for production is listed in 
Table E] and Ell and plotted in Fig. [HI The A{4>) data 
were fitted with a sine curve, Eq. (|lll) . to obtain An with 
4>a = 0. The mean if values were determined according 
to section HV Al Statistical uncertainties are shown as er- 
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TABLE V: The results of the xf dependence of An for neu- 
tron production in the ZDC trigger sample of p+p collisions 
at -^5=200 GeV. First and second uncertainties show statis- 
tical and pT-correlated systematic uncertainties, respectively. 
Scale uncertainties from the asymmetry measurements and 
the beam polarization are not included. 



X 2 /ndf 



-0.776 

-0.682 

-0.568 

0.568 

0.682 

0.776 



-0.0059±0.0252±0.0095 
-0.0219±0.0255±0.0035 
-0.0050±0.0303±0.0076 
-0.0503±0.0263±0.0076 
-0.0625±0.0221±0.0035 
-0.0772±0.0217±0.0095 



11.6/7 
6.833/7 
9.252/7 
7.012/7 
2.68/7 
5.38/7 



TABLE VI: The results of the xf dependence of An for 
neutron production in the ZDC®BBC trigger sample of p+p 
collision at -^5=200 GeV. First and second uncertainties show 
statistical and pT-correlated systematic uncertainties, respec- 
tively. Scale uncertainties from the asymmetry measurements 
and the beam polarization are not included. 



X 2 /ndf 



-0.749 

-0.664 

-0.547 

0.547 

0.664 

0.749 



0.0035±0.0117±0.0082 
-0.0093±0.0106±0.0037 
-0.0033±0.0115±0.0096 
-0.0629±0.0097±0.0096 
-0.0657±0.0090±0.0037 
-0.0667±0.0099±0.0082 



2.672/7 
2.915/7 
6.783/7 
13.27/7 
5.425/7 
5.003/7 



ror bars and ^-correlated systematic uncertainties are 
shown as brackets. Scale uncertainties from the asym- 
metry measurements and the beam polarization are not 
included. 

We observe significant negative An for neutron pro- 
duction in the positive xp region and with no energy 
dependence within the uncertainties, both for inclusive 
neutron production and for production including a beam- 
beam interaction requirement. No significant backward 
neutron production asymmetry is observed. 



V. DISCUSSION 

The measurement of the cross section for the p+p pro- 
duction of neutrons at y/s—200 GeV has been presented 
here, and it is consistent with xf scaling when compared 
to ISR results. These cross sections are described by the 
OPE model in Regge calculus pl-Ucj. Therefore, the ob- 
served large asymmetry for neutron production at RHIC, 
as presented in [l[ and here, may arise from the inter- 
ference between a spin-flip amplitude due to the pion 
exchange and nonflip amplitudes from other Reggeon ex- 
changes. So far our knowledge of Reggeon exchange com- 
ponents for neutron production is limited to the pion. 
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FIG. 17: The xf dependence of An for neutron produc- 
tion in the ZDC trigger sample (top) and for the ZDC®BBC 
trigger sample (bottom). The error bars show statistical un- 
certainties and brackets show pT-correlated systematic uncer- 
tainties. Systematic scale uncertainties listed in Table IWl and 
polarization scale uncertainties are not included. 



Under the OPE interpretation, the asymmetry has sen- 
sitivity to the contribution of all spin nonflip Reggeon 
exchanges, even if the amplitudes are small. Recently 
Kopeliovich et al. calculated the An of forward neutron 
production from the interference of pion and Reggeon ex- 
changes, and the results were in good agreement with the 
PHENIX data 

We can also discuss our results based on the meson- 
cloud model [Hj]. This model gives a good description 
for the result from a Drell-Yan experiment at FNAL, 
E866[23|]. In this model, the Drell-Yan process is gen- 
erated by the interaction between the d quark in one 
proton and the d quark in the 7r + of p — > mr + state 
for other proton. In this model the neutron should be 
generated with very forward kinematics, possibly simi- 
lar to the kinematics of the results presented here. The 
meson-cloud model was successfully applied to neutron 
production in the ISR experiment [8j and we expect it is 
applicable to our An and cross section measurements for 
higher energy p+p collisions. 
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VI. CONCLUSION 

We have measured the cross section and single trans- 
verse spin asymmetry, An, for very forward neutron 
production in polarized p+p collisions at -^=200 GeV. 
The results from the PHENIX experiment at RHIC were 
based on a zero degree hadronic calorimeter (ZDC) aug- 
mented by a shower maximum detector, covering neu- 
tron production angles to 6>„=2.2 mrad. A large An for 
neutron production had been observed in a polarimeter 
development experiment at RHIC, using an electromag- 
netic calorimeter to identify neutrons, with coarse neu- 
tron energy resolution [l]]. The PHENIX experiment then 
outfitted existing ZDC detectors to act as polarimctcrs 
to monitor the beam polarizations and polarization direc- 
tions at the experiment. The results presented here are 
based on studies with the ZDC polarimeter, which due to 
a much better measurement of the neutron energy, pro- 
vide first measurements of the neutron production cross 
section at RHIC energy, and the dependence of An on 
the neutron energy. 

The measured cross section is consistent with x p scal- 
ing from ISR results. Within uncertainties, the observed 
An were consistent with the previous result at RHIC [1] 
and for xf > 0.45 (the region measured by this exper- 
iment) no significant xf dependence was observed. We 
also present measured An for neutrons produced back- 
ward from the polarized beam. These results are consis- 
tent with zero. 

The cross sections for large xp neutron production, as 
well as those in e + p collisions at HERA, are largely re- 
produced by a one pion exchange model (OPE). Using 
this model, the observed large asymmetry for the neu- 
tron production would be considered to come from the 
interference between a spin-flip amplitude due to the pion 
exchange and nonflip amplitudes from other Rcggeon ex- 
changes. On the basis of the OPE model, the large neu- 
tron An would have sensitivity to the contribution of 
other Reggeon exchanges. 

Future measurements of neutron production cross sec- 
tions and asymmetries will include analysis of RHIC runs 
at ^=62.4 GeV and at 500 GeV. The measurements at 
different center of mass energies will probe the x f and pt 
dependence for neutron production at fixed, very forward 
production angles 9 n <2.2 mrad. 
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Appendix A: The study of the beam axis on the 
detector geometry 

The ZDC center was aligned to the beam axis at the 
beginning of the 2003 run. We assumed that the beam 
axis was on the ZDC center in this analysis of 2005 data. 
We used two approaches to estimate the beam and ZDC 
alignment. Peripheral neutrons from a heavy ion run 
just prior to the polarized proton run gave centers of 
x = 0.28 ± 0.01 cm and y = -0.07 ± 0.01 cm at the 
south ZDC. The center of the asymmetry An was also 
used to determine the center of the beam axis at the 
ZDC, since An must be zero at zero production angle. 
We used the ZDC<g)BBC trigger sample in this analysis. 
The asymmetry was measured for a vertically polarized 
beam to obtain the center in x, and for a special run 
with horizontally polarized beam to obtain the center 
in y. The results were x = +0.46 ± 0.08 cm and y — 
— 1.10±0.14 cm. The results of the two techniques agreed 
reasonably for x, and did not agree for y. 

The beam axis shifts that we observed were consid- 
ered as systematic uncertainties for the results. The un- 
certainties were determined from variations of the cross 
section and asymmetry obtained by moving the center of 
acceptance while keeping the same cut region (for exam- 
ple, r<2 cm for the cross section analysis). 
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Appendix B: Energy Unfolding 

The measured neutron energy with the ZDC is smeared 
by the energy resolution. For the extraction of the origi- 
nal energy distribution, it is necessary to unfold the mea- 
sured energy distribution. We use an unfolding method 
proposed in 

We assume that the initial distribution x(E) is smeared 
to the measured distribution y(E') and this smearing is 
described by a linear combination. Their relation can be 
given by a transition matrix A(E' , E) as, 



y(E')=A(E',E)x(E). 



(Bl) 



or y = Ax. 



If the smearing effect is large, the result is very sensi- 
tive to a small change of A. It can be discussed using an 
orthogonal decomposition. The matrix A is diagonalized 
into D with a transformation matrix U, 



'= Db, 



(B2) 



where D = U~ 1 AU, and c = U y and b = U x are 
new vectors transformed from y and x, respectively. The 
diagonal elements of the matrix D are the eigenvalues 
Xj of the matrix A. Each of the coefficients bj and Cj 
in c = Db is transformed independently of any other 
coefficient by using eigenvalue A,,-, 



x10 b 



Cj = Xj- bj. 



(B3) 
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FIG. 18: Energy distributions before and after the energy 
unfolding. The unfolding was performed so that statistics 
was conserved. 



In order to perform the unfolding, the coefficients Cj 
have been affected by statistical fluctuations of the ele- 
ments of measured vector y. The bj which includes the 
information of initial vector x is obtained by bj = Cj/Xj. 
The statistical fluctuation of the Cj amplified in the case 
of small eigenvalue Xj , resulting in instability. Reason- 
able result can be obtained by cutting the Cj which has 
a large statistical uncertainty. 

First, the coefficients Cj were calculated. Three sets 
of the transition matrix A, which have the same energy 
resolution but different initial energy distributions, were 
prepared with a simulation to check the statistical er- 
ror propagation of the Cj . Initial shapes were prepared 
to increase, be flat and decrease as a function of xp. 
These shapes are close to the cross sections at px ~ 0.0 
GeV/c, 0.2 GeV/c and 0.4 GeV/c in the ISR results. En- 
ergy spectra before and after the unfolding are plotted in 
Fig- UHl Horizontal axis is changed to Xp by Eq. |(5J). 
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